Monte Carlo simulation of equilibrium Ll ordering in FePt nanoparticles 
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First, second and third nearest neighbor mixing potentials for FePt alloys, were calculated from 
first principles using a Connolly- Williams approach. Using the mixing potentials obtained in this 
manner, the dependency of equilibrium Llo ordering on temperature was studied for bulk and for a 
spherical nanoparticle with 3.5nm diameter at equiatomic composition by use of Monte Carlo sim- 
ulation and the analytical ring approximation. The calculated order-disorder temperature for bulk 
(1495-1514 K) was in relatively good agreement (4% error) with the experimental value (1572K). 
For nanoparticles of finite size, the (long range) order parameter changed continuously from unity 
to zero with increasing temperature. Rather than a discontinuity indicative of a phase transition we 
obtained an inflection point in the order as a function of temperature. This inflection point occurred 
at a temperature below the bulk phase transition temperature and which decreased as the particle 
size decreased. Our calculations predict that 3.5nm diameter particles in configurational equilibrium 
at 600° C (a typical annealing temperature for promoting Llo ordering) have an Llo order parameter 
of 0.83 (compared to a maximum possible value equal to unity). According to our investigations, 
the experimental absence of (relatively) high Llo order in 3.5nm diameter nanoparticles annealed 
at 600° C or below is primarily a problem of kinetics rather than equilibrium. 



I. INTRODUCTION 

Self-assembled, monodispersed FePt nanoparticles are 
being intensively investigated for possible future applica- 
tion as an ultra-high density magnetic storage medium. 
In order to be useful as a storage medium, these parti- 
cles, because of their extremely small volume, V, must 
have sufficiently high magnetic anisotropy, K Ul to with- 
stand thermal fluctuations of the direction of magnetiza- 
tion. This requires values of the thermal stability factor, 
(K u V)/(k%T), of approximately 50. The particles are 
usually produced by a "hot soap" process that yields a 
disordered fee solid solution alloy (e.g. Ref. 0). Such 
particles are not useful for information storage because 
they are superparamagnetic at room temperature due to 
their low magnetic anisotropy. 

Typically, the particles are annealed at a tempera- 
ture T ~ 600° C in order to induce an ordered Ll 
phase 2 * 3 -. The layered Llo phase* is known from stud- 
ies of bulk alloys to have an extremely high magnetic 
anisotropy (K u = 7 x 10 7 erg/cm 3 ). This value of mag- 
netic anisotropy would provide a sufficiently large ther- 
mal stability factor to make 3.5nm diameter particles vi- 
able for information storage. 

Unfortunately, it appears to be difficult to a achieve a 
high degree of long range atomic order in FePt nanopar- 
ticles with <4nm diameter by annealing at T <600°C 
(e.g. Ref 2). One can consider two possible reasons for 
the fact that it has not been possible to obtain well or- 
dered small particles. First, the observed order may be 
low because the particle is not in its equilibrium state due 
to the slow kinetics at low temperatures. Alternatively, 
the equilibrium order itself may be low even at relatively 
low temperatures because of the small size of nanopar- 



ticles. The latter explanation was suggested in Ref. 2. 
There, the order-disorder phase transition temperature 
was estimated to decrease with decrease of particle size. 
For particle sizes less than 1.5 nm in diameter, the phase 
transition temperature was found to be below the typical 
annealing T ~ 600°C. Therefore, particles of diameter 
less than 1.5 nm were predicted to have no long range 
order in their equilibrium state at 600°C. This explana- 
tion is in qualitative agreement with experiment. The 
difference between the experimental (4mn) and theoret- 
ical (1.5nm) critical size for disappearance of Llo order 
at 600°C was attributed to the neglect of nanoparticle 
surface effects. 

From our point of view, however, the results obtained 
in Ref. :2j require verification because of the limitations of 
the theoretical models used in that study. Namely, the 
interatomic potentials in alloys usually are much more 
complicated and long-ranged than the nearest neigh- 
bor Lennard- Jones model that was used. In addition, 
the order-disorder phase transition temperature was es- 
timated in Ref. y by comparing the free energies of 
completely ordered and completely disordered states; 
whereas in reality, the ordered state approaches (with 
increasing temperature) the phase transition point be- 
ing not completely ordered. Also, the disordered state 
would be expected to approach the phase transition (with 
decreasing temperature), not with a completely random 
atomic distribution but with an atomic distribution that 
has substantial short range order. Moreover, it is known 5 
that there is no formal phase transition in a finite system. 

In the present paper we utilize first principles cal- 
culations (VASP code&) together with the Connolly- 
Williams 7 " method and Monte Carlo (MC) simulations 
(utilizing the Metropolis algorithm 8 ) to study the tem- 
perature dependence of equilibrium L1q order in a 
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spherical FePt nanoparticle with 3.5nm diameter and 
equatomic composition (c = 0.5). 



II. RESULTS 

We consider an Fe-Pt alloy in the framework of the 
commonly used two-component lattice gas model. In 
such a model 9 , two types of atoms are distributed over 
the sites of a rigid crystal lattice. The atoms are al- 
lowed to be situated only at the crystal lattice sites and 
each site can be occupied by only one atom. The atoms 
interact through the lattice potentials (so-called mixing 
potentials) and can exchange their positions according to 
Gibbs statistics. 

We used the Connolly- Williams 7 method to calculate 
the mixing potentials. Within this method, the ener- 
gies of several periodic atomic distributions (i.e. long- 
range ordered structures called superstructures; for ex- 
ample Llo) are calculated by first principles methods. 
Then the mixing potentials are determined by the best 
fit to those energies. We considered twenty three lin- 
early independent Fe-Pt superstructures of the same 
equiatomic composition c = 0.5. First principles calcu- 
lations were performed within the local-density approxi- 
mation to density- functional theory, using the VASP pro- 
gram package. 6 ' All superstructures were totally relaxed 
including shape and volume relaxation of the unit cell 
and individual displacements of atoms within the unit 
cell. An 8x8x8 mesh of fc-points in the full Brillouin 
zone was employed. 

The Llo superstructure was included in our first set of 
first principles calculations. In this case, after atom po- 
sition relaxation, we obtained 3.848A and 3.77lA for a 
and c lattice parameters of the corresponding tetragonal 
lattice, respectively. For comparison the experimental 
values are 3.847A and 3.715A4a In addition, our calcu- 
lated results showed the Llo ferromagnetic superstruc- 
ture to be more stable (i.e. has lower energy) than the an- 
tiferromagnetic one in accordance with experiment. We 
believe that this good correspondence between theoret- 
ical and experimental results confirms the adequacy of 
our VASP first principles calculations. 

By applying the Connolly- Williams method, we ob- 
tained 0.08769 eV, -0.03946 eV and 0.01585 eV for the 
first, second and third nearest neighbor pair mixing po- 
tentials, respectively. The average accuracy with which 
we fit the energy of the twenty-three superstructures 
within the Connolly- Williams method was 1.14% per one 
structure. 

To verify the calculated values of mixing poten- 
tial, we calculated the phase transition temperature in 
the bulk FePt alloy using these values. As a result 
we obtained 1495 K and 1514 K within the analyti- 
cal ring approximation^ and MC simulation, respec- 
tively. The close correspondence of these values to the 
experimental^ one of 1572 K (4% error), demonstrates 
the adequacy of the calculated mixing potential. 
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FIG. 1: The temperature dependence of the FePt equilib- 
rium Llo order parameter r\ in the cases of bulk ("bulk") 
and a spherical nanoparticle with 3.5nm diameter ("sphere") 
at equiatomic composition c = 0.5. Two results for bulk 
were obtained by Monte Carlo simulation ("MC") for par- 
allelepiped sample containing N = 216000 atoms and within 
the analytical ring approximation— ("ring"). At simulation, 
the starting configuration for each temperature was chosen 
to be the completely ordered one. We applied free and pe- 
riodic boundary conditions in cases of spherical nanoparticle 
and parallelepiped, respectively. For the case of a nanoparti- 
cle at 378°C, 528°C and 600°C, the error bars correspond to 
dispersion of r/ due to the thermodynamic fluctuations. 

To investigate long range order in spherical nanoparti- 
cles, we used the calculated mixing potentials in MC sim- 
ulations to determine the temperature dependence of the 
equilibrium Llo order parameter rj in the case of spherical 
FePt nanoparticlcs with 3.5nm diameter and equiatomic 
composition c = 0.5. The results are presented in Fig. ^ 

We define the equilibrium Llo order parameter rj as the 
statistical average of the maximum value among three ab- 
solute values of "directional" order parameters rj Xl rj yi rj z : 

r} = (max{ |r? x |,|%|, \rj z | }) MC , (1) 

where r/i (i — x, y, z) is defined as difference between 
the Fe atom concentrations at odd and even crystal 
planes perpendicular to i-th direction, (. . .} M q is the sta- 
tistical average over the MC steps. We chose such a defi- 
nition of rj because of the equivalence by symmetry of the 
x, y, and z directions of Llo order. In addition, one can 
obtain an equivalent structure (at c = 0.5) by changing 
the sign of rji, which results in the exchange of Fe and 
Pt atoms producing a configuration that is equivalent by 
symmetry to the original one. During MC simulation, 
we observed fluctuations that cause transformations be- 
tween these equivalent states (i.e. fluctuations in the sign 
and direction of 77)— This is in addition to the usual sta- 
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tistical fluctuations within one such state. The Llo order 
parameter 77, defined in Eq. (1} takes into account the 
fluctuation induced transformations between the equiva- 
lent states^. 



III. CONCLUSIONS AND DISCUSSION 

From Fig. n° ne may conclude the following. The ring 
approximation (which corresponds to bulk, i.e. an infi- 
nite sample) clearly shows a phase transition when the 
order parameter r\ drops to zero. Strictly speaking, in 
both of the cases considered here of finite size samples 
(sphere and parallelepiped) there is no phase transition 
in accordance with a general theorem^. The order pa- 
rameter r\ continuously changes from unity to zero with 
increasing temperature and instead of a phase transition 
we obtain an inflection point in the rj(T) curve. In the 
case of the parallelepiped with 216000 atoms, the inflec- 
tion point is very similar to the phase transition^ 

Our calculations predict that 3.5nm diameter particles 
in configurational equilibrium at 600° C would have an 
order parameter r\ = 0.83 (compared to a maximum pos- 
sible value of unity). Therefore, annealing at 600°C will 
not yield perfect order for 3.5nm diameter particles. Ap- 
proximately 17% of the atoms will be on the wrong sub- 
lattices, even in equilibrium. The dispersion of 77 due to 
the thermodynamic fluctuations is comparatively small 
(2.5%) near annealing T = 600°C. 

According to our investigations, the experimental ab- 
sence of (relatively) high order in nanoparticles below 
600° C is primarily a kinetic problem rather than an equi- 
librium one. It should be noted that to rapidly obtain the 
correct equilibrium state, we used simplified kinetics in 
our MC simulation^. Namely, we allowed any two ran- 
domly chosen atoms to exchange their positions without 
an additional diffusion barrier. In a real alloy, the main 



mechanism of atomic diffusion is much slower because 
it consists in exchange the positions between atoms and 
their nearest neighbor vacancies through energy barriers. 
Moreover, at each temperature we started the simulation 
from the completely ordered state, whereas the actual 
nanoparticles are initially prepared in disordered state 
and transformation from the disordered to the ordered 
state may be much slower than the reverse one, espe- 
cially at low temperatures. Nevertheless, even with our 
simplified kinetics, we observed a slowing down problem 
in approaching the equilibrium ordered state at low tem- 
peratures. In real nanoparticles this problem must be 
much worse. Kinetic acceleration methods such as irra- 
diation and/or addition of other types of atomsi^ may be 
useful in accelerating the formation of long range order. 

In our study we used mixing potentials obtained for 
infinite bulk alloys and used free boundary conditions to 
simulate the equilibrium configuration of finite size par- 
ticles. The presence of the surface will change the atomic 
potentials in the near-surface region in comparison with 
bulk potentials. Analytical estimation of such surface ef- 
fects is not straightforward and will be done elsewhere 14 . 
In reality, the problem of the effect of the surface on 
the interatomic exchange potentials is even more compli- 
cated because the nanoparticles of most current interest 
are likely to have unknown atoms and molecules attached 
to their surfaces. 
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